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' fitfinder Macro 

1 Macro recorded 9/29/00 by Michael J. Corey 

' Keyboard Shortcut: Ctrl+Shift+A 

Rem initialize rate and correlation arrays 
Rem 

Rem Fit (rate) array, allowing for future expansion 
Dim fit_ar(5000, 96) 

Rem Correlation array, which is searched for the best correlation later 
Dim cor_ar(1000, 96) 
tempi = 1 

Rem get control data from EGG file 
tp_no = Range(" A4") 
sample_no = Range("A3") 
Rem This initializes the time array 
For ix = 61 To 60 + tp no 
Cells(ix, 1). Value - Cells(ix - 48, 1) 
Next ix 
Rem 

Rem This begins the outer loop. Each "type" of fit (e.g., 

Rem the first three points, or points 7-12) is done for all 

Rem samples and the results are stuffed into the rate and 

Rem correlation arrays. Then we go back and do the next "type" 

Rem of fit. 

Rem 

temp2 = 0 

For stix - 13 To (tp_no + 12) 

Rem Initialize the number of points to fit to 3, since there 
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Rem is little point in getting a linear fit of 2 points. 

pts_to_fit = 4 

Rem 

Rem This inner loop controls the endpoint of the fit. 
Rem 

For eix = (stix + 3) To (tp_no + 12) 
Rem 

Rem This loop repeats the fit for each data- set (i.e., sample). 
Rem 

temp2 = temp2 + 1 
For sample ix = 2 To sample no + 1 
o Rem This copies the data from the raw data matrix to a fixed position. 

M- Range( ,f B60"). Select 

m For pts_jx = 1 To pts_to_fit 

Q ActiveCell.Offset(l, 0). Select 

o ActiveCell. Value = Cells(stix + pts_Jx - 1 ? sample_ix) 

fi: Nextpts ix 

: =ff Rem The row index into the arrays is calculated beforehand 
fij Rem to make coding easier. 
Rem 

Rem Case routine on number of points to fit 
Rem 

If pts_to_fit = 4 Then GoTo fit4 
If pts to fit = 5 Then GoTo fit5 
If pts_to_fit = 6 Then GoTo fit6 
If pts_to_fit = 7 Then GoTo fit7 
If ptsjojit = 8 Then GoTo fit8 
If ptsjo fit = 9 Then GoTo fit9 
If pts_to_fit = 10 Then GoTo fit 10 
fit4: 
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Range("B2").FormulaRlCl = "=LinEst(r61c2:r64c2,r61cl:r64cl)" 
Range("C2").FormulaRlCl = "=correl(r61c2:r64c2,r61cl:r64cl)" 
GoTo fitend 
fit5: 

Range( H B2").FormulaRlCl = "=LinEst(r61c2:r65c2,r61cl:r65cl)" 
Range("C2").FormulaRlCl = M =correl(r61c2:r65c2,r61cl:r65cl)" 
GoTo fitend 
fit6: 

Range("B2").FormulaRlCl = "=LinEst(r61c2:r66c2,r61cl:r66cl)" 
Range("C2 H ).FormulaRlCl = "=correl(r61c2:r66c2,r61cl:r66cl)" 
GoTo fitend 
fit7: 

Range("B2").FormulaRlCl = "=LinEst(r61c2:r67c2,r61cl:r67cl)" 
Range("C2").FormulaRlCl = "=correl(r61c2:r67c2,r61cl:r67cl)" 
GoTo fitend 
fit8: 

Range("B2").FormulaRlCl = "=LinEst(r61c2:r68c2,r61cl:r68cl)" 
Range("C2").FormulaRlCl = "=correl(r61c2;r68c2,r61cl:r68cl)" 
GoTo fitend 
fit9: 

Range("B2").FormulaRlCl = "=LinEst(r61c2:r69c2,r61cl:r69cl) M 
Range("C2 ,, ).FormulaRlCl = "=001X610-6102x6902x6101x6901)" 
GoTo fitend 
fitlO: 

Range("B2").FormulaRlCl = "=LinEst(r61c2:r70c2,r61cl:r70cl)" 
Range("C2").FormulaRlCl = "=correl(r61c2:r70c2,r61cl:r70cl)" 
GoTo fitend 
fitend: 

fit_ar(templ, sampleix - 1) = Range("B2") 

cor_ar(templ, sample_ix - 1) = Range("C2") 

Rem The next four lines display interim data in various parts of 



Rem the spreadsheet, mostly to make debugging easier. 
temp3 = Range("B2") 
temp4 = Range("C2") 

Cells(temp2 + 30, sample Jx) = Range("B2").Value 

Cells(temp2 + 90, sample Jx) = Range("C2"). Value 

Next sample_ix 

Rem increment the fit counter 

tempi = tempi + 1 

Rem increment the length of range to fit 

pts to fit = pts Jo Jit + 1 

Rem current limit is ten points in the fit range 

If pts Jo Jit = 1 1 Then GoTo endthiseix 

dumberlabel: 

Next eix 

Next stix 

GoTo skipend 

endthiseix: 

eix = tp_no + 13 

GoTo dumberlabel 
Rem 

Rem Now the routine that determines the best fit. 

Rem I will scan from the end of the array so that if there 

Rem is a tie, it will pick up the run with the most data points. 

Rem 

skipend: 

Cells(8, l).Value - "max correl" 
Cells(9, l).Value - "rate" 
Cells(10, l).Value = "notes" 
For samplejx = 1 To sample_no 
no_of Jits = (tp_no - 2) * (tp_no - 3) / 2 

Rem Note: the previous calculation assumes the fit begins with 4 points. 
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Rem In the current code the maximum size of fit is 10 points, and this 
Rem formula doesn't take that into account-it would have to be changed 
Rem to accommodate fits of more than 10 points, but the principle would 
Rem be the same, 
equalflag = 0 

Dim correlmax As Variant 
Dim tempcorr As Variant 
correlmax = -1.1 
fitmax = -1 

For fitix = no_of_fits To 1 Step -1 
tempcorr = cor_ar(fitix, sample_ix) 

wis: 

□ Rem Note: the correl function can yield a DIV/0 error 

H i; 

sj Rem when the correlation is somehow strange, for example 
H ; Rem in negative-control reactions. 

E : i: 

111 Rem The next statement skips the cases in which the error 

7" Rem appears. 

O If IsNumeric(tempcorr) Then 

Q GoTo nexterl 

p( Else: GoTo nexter 

fU End If 

nexterl : 

If tempcorr = correlmax Then equalflag = 1 

If tempcorr > correlmax Then GoTo prenexter 

GoTo nexter 

prenexter: 

equalflag = 0 

fitmax = fitix 

correlmax = tempcorr 

Rem This next complicated routine deconvolves the triangle number representing 
Rem the mapping between the index into the fit array and the actual time range that 
Rem was used—for example, using 9 timepoints, a fitix of 3 would be the third fit, 
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Rem which would be timepoints 1 to 6; a fitix of 10 would be timepoints 2 to 8. 

Rem obviously this routine can be made smaller by having some kind of "for" loop 

Rem instead of all the if statements. That would also allow it to handle more than 

Rem ten timepoints, which is currently the limit of the whole program. 

If fitix > tp_no - 3 Then GoTo tryx2 

fitlo = 1 

fithi = fitix + 3 

GoTo nexter 

tryx2: 

If fitix > 2 * tp_no - 7 Then GoTo tryx3 
fitlo = 2 

f? fithi = fitix + 7 - tp_no 
O GoTo nexter 

""t-S 

li tryx3: 

W If fitix > 3 *tp_no- 12 Then GoTo tryx4 

p fitlo = 3 

p fithi = fitix + 12 - tpjio * 2 
W GoTo nexter 
ffl tryx4: 

H If fitix > 4 * tp_no - 1 8 Then GoTo tryx5 

fitlo = 4 

fithi = fitix + 18 -tp_no* 3 

GoTo nexter 

tryx5: 

If fitix > 5 * tp_no - 25 Then GoTo tryx6 
fitlo = 5 

fithi = fitix + 25 - tp_no * 4 

GoTo nexter 

tryx6: 

If fitix > 6 * tp_no - 33 Then GoTo tryx7 
fitlo = 6 
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fithi = fitix + 33 - tp_no * 5 
GoTo nexter 

tryx7: 

If fitix > 7 * tp_no - 42 Then GoTo tryx8 
fitlo = 7 

fithi = fitix + 42 - tp_no * 6 

GoTo nexter 

tryx8: 

If fitix > 8 * tp_no - 52 Then GoTo tryx9 
fitlo = 8 

fithi = fitix + 52 - tp_no * 7 

GoTo nexter 

tryx9: 

If fitix > 9 * tp_no - 63 Then Stop 

fitlo = trno - 1 1 

fithi = tp_no- (fitix -37) 

GoTo nexter 

nexter: 

Next fitix 

Rem Display the hard-won best correlation, associated rate, and time range used. 

Cells(8, sample_ix + 1) = correlmax 

Cells(9, sample_ix + 1) = fit_ar(fitmax, sample_ix) 

Cells(10, sample_ix + 1) = "fit " + CStr(fitlo) + " to " + CStr(fithi) 

Next sample_ix 

End Sub 
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1 

9 
9 

MJC5.RF. 20.11.2000 
2.01 1.24 



14:15:06 
3 



Time A3 A4 A5 B3 B4 B5 C3 C4 C5 



0 


142337 


198427 


282680 


10882 


17286 


50090 


120279 


137692 


136200 


31.2 


154969 


222070 


338398 


11467 


18024 


51071 


119729 


136967 


135436 


61.3 


166526 


244896 


388966 


11583 


18137 


51124 


118809 


135772 


134596 


91.3 


177288 


266025 


435282 


11657 


18217 


51120 


117688 


134883 


133384 


121.3 


187991 


286313 


478786 


11824 


18186 


50952 


116941 


133818 


132023 


151.3 


198051 


306084 


519934 


11978 


18296 


50721 


115978 


132570 


131284 


181.3 


208351 


325526 


557892 


12006 


18340 


50636 


115283 


131989 


130199 


211.3 


218437 


343973 


593623 


12062 


18326 


50507 


114468 


130781 


129117 


241.3 


228325 


362271 


627305 


12259 


18279 


50344 


114006 


130131 


128443 



2413 AppthJ-tX 3 

1 -31.614891 -0.995236 ' ' 

9 
9 

MJC5.RPT 20.11.2000 14:15:06 
2.01 1.24 3 



max correl 0.9999951 0.999997 0.99984177 0.9902173 0.956535 0.8042586 -0.98781 -0.98966 -0.99095 
rate 331.93168 652.21619 1435.7657 4.1695772 1.919422 10.411267 -28.554 -28.023 -34.3462 

notes fit 6 to 9 fit 4 to 7 fit 3 to 6 fit 2 to 6 fit 2 to 7 fit 1 to 4 fit 1 to 4 fit 6 to 9 fit 1 to 5 



Time A3 A4 AS B3 B4 B5 C3 C4 C5 



0 


142337 


198427 


282680 


10882 


17286 


50090 


120279 


137692 


136200 


31.2 


154969 


222070 


338398 


11467 


18024 


51071 


119729 


136967 


135436 


61.3 


166526 


244896 


388966 


11583 


18137 


51124 


118809 


135772 


134596 


91.3 


177288 


266025 


435282 


11657 


18217 


51120 


117688 


134883 


133384 


121.3 


187991 


286313 


478786 


11824 


18186 


50952 


116941 


133818 


132023 


151.3 


198051 


306084 


519934 


11978 


18296 


50721 


115978 


132570 


131284 


181.3 


208351 


325526 


557892 


12006 


18340 


50636 


115283 


131989 


130199 


211.3 


218437 


343973 


593623 


12062 


18326 


50507 


114468 


130781 


129117 


241.3 


228325 


362271 


627305 


12259 


18279 


50344 


114006 


130131 


128443 



Li; 
Jj I 



n 



383.01808 
375.45678 
367.59442 
362.18593 
357.769 

353.8321 
361.30767 
355.60397 
352.21878 

349.2975 
346.44278 
346.31774 
344.94916 
343.31575 
341.26961 
339.64045 
339.14603 

337.6527 
334.29926 
333.82595 
331.93168 



742.25889 
726.00556 
711.83165 
699.95175 
687.9326 
677.13031 
703.59187 
692.02935 
682.46681 
671.84166 
662.26498 
670.59796 
665.11777 
656.42933 
648.35926 
652.21619 
644.59517 
637.48812 
632.98969 
627,07719 
615.18674 



1672.7546 
1616.26211 
1565.31615 
1515.08105 

1467.2062 
1421.93685 
1538.10089 
1496.52598 
1451.78988 
1407.96718 

1365.9922 

1435.7657 
1396.06535 
1355.58138 
1316.25617 
1345.55402 
1307.80713 

1270.7113 
1258.37217 
1225.00586 

1177.2935 



8.066686 
6.8755017 
6.2845871 
5.5029586 
4.8744994 
4.6925229 

3.764263 
4.1695772 
3.8289985 
3.4978006 
3.6133597 
4.4406335 
3.8544857 
3.4037896 
3.5779737 
3.9581084 
3.2838643 
3.5532585 
2.4486011 
3.1524772 
2.9472259 



9.607484 
6.621386 
5.339885 
4.466975 

3.69283 

2.98143 
1.870181 
1.962817 
1.919422 
1.650081 

1.27085 
1.467754 
1.601985 
1.348177 
0.946909 

1.56841 
1.226795 

0.73285 
1.534317 
0.721154 
-0.20762 



10.411267 
5.9133514 
2.6869172 
0.9391404 
-0.226583 
-1.117821 
-1.174205 
-2.866657 
-3.351375 
-3.626635 
-3.891538 
-4.513374 

-4.53421 
-4.507104 
-4.595845 

-5.53882 
-5.097408 
-5.017277 
-4.679316 
-4.728085 
-4.139197 



-28.554 
-28.7757 
-29.3101 
-28.8491 
-28.5167 
-27.4908 
-31.2029 
-30.9569 
-29.7838 
-29.0835 
-27.6568 
-30.4087 
-28.9595 
-28.2868 
-26.6744 
-26.8971 
-26.7511 
-25.1032 
-26.7022 
-24.3964 
-22.1495 



-31.6295 
-32.4664 
-34.0004 
-32.9962 

-33.178 
-32.4995 
-34.0106 
-35.5059 
-33.6553 
-33.6925 
-32.7119 
-35.0756 
-32.6337 
-33.0073 
-31.9358 
-32.6838 
-33.1527 
-31.6635 
-31.8953 
-30.2849 

-28.023 



-30.524 
-34.3462 
-34.0775 

-34.203 
-34.41 
-33.7145 
-37.6282 
-35.9209 
-35.4619 
-35.3611 
-34.2507 
-37.1803 
-36.0042 
-35.6729 
-34.1494 
-33.8872 
-34.2299 
-32.6651 
-32.2178 
-30.8044 
-31.6149 



131284 
130199 
129117 
128443 
128443 
128443 
128443 
128443 
128443 



0.999633 
0.9996012 
0.9994528 
0.9994593 
0.9994765 
0.9994797 
0.9999576 
0.9998471 
0.9998482 
0.9998444 

0.999826 
0.9999425 
0.9999627 
0.9999631 
0.9999464 
0.9999861 
0.9999919 
0.9999821 
0.9999258 
0.9999614 
0.9999951 



0.9998448 
0.9996634 
0.9995317 
0.9994544 
0.9993158 
0.999207 
0.999843 
0.9997762 
0.9997355 
0.9996188 
0.9995276 
0.9999801 
0.9999553 
0.9998555 
0.999777 
0.999997 
0.9999274 
0.9998735 
0.9999557 
0.999932 
0.9999875 



0.99947062 
0.99909356 
0.99873209 
0.99823055 
0.99767264 
0.99708114 
0.99968876 
0.99944197 
0.99902011 
0.99853484 
0.99801335 
0.99984177 
0.99950457 
0.99909758 
0.99865159 
0.99975498 
0.99944582 
0.99908896 
0.99973904 
0.99948407 
0.99983614 



0.8986993 
0.9144503 
0.9339036 
0.9309838 
0.9267276 
0.9416045 
0.9877688 
0.9902173 
0.9881171 
0.9837782 

0.988755 
0.9869261 

0.979931 
0.9739874 
0.98309 
0.9659756 
0.9554568 

0.973789 
0.9436501 
0.9586433 
0.9126798 



0.882256 
0.813995 
0.809147 
0.806686 
0.787346 
0.746819 
0.866925 
0.931044 
0.956535 

0.94338 
0.870883 
0.863685 

0.93197 
0.924294 
0.811485 
0.869944 
0.862655 
0.679448 
0.862855 
0.570086 
-0.29348 



0.8042586 
0.6398954 
0.37913 
0.1635285 
-0.046029 

-0.24955 
-0.574589 

-0.80943 
-0.894652 
-0.935272 
-0.956765 
-0,932756 
-0.964636 

-0.97896 
-0.986851 
-0.987001 
-0.988493 
-0.992962 
-0.980615 

-0.99023 
-0.989244 



-0.98781 
-0.99386 
-0.99637 
-0.99743 
-0.99811 
-0.99695 
-0.99746 
-0.99866 
-0.99813 
-0.99828 
-0.99607 
-0.99774 
-0.99748 

-0.9981 
-0.99561 
-0.99804 
-0.99898 
-0.99645 

-0.9985 
-0.99459 
-0.99499 



-0.9954 
-0.99746 

-0.9973 
-0.99743 
-0.99827 
-0.99822 
-0.99888 
-0.99858 
-0.99706 
-0.99816 
-0.99779 
-0.99648 
-0.99512 
-0.99717 
-0.99716 
-0.99116 
-0.99554 
-0.99571 
-0.99164 
-0.99388 
-0.98966 



-0.99285 
-0.99095 
-0.99465 
-0.99665 
-0.99775 
-0.99786 
-0.99368 
-0.99535 
-0.99714 
-0.99819 
-0.99766 
-0.99391 
-0.99617 
-0.9977 
-0.99676 
-0.99463 
-0,99728 
-0.99675 
-0.99548 
-0.99644 
-0.99524 
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